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Abstract 

The production rate of heavy Majorana neutrinos is relevant for models of thermal 
leptogenesis in the early Universe. In the high temperature limit the production can 
proceed via the 1 -H- 2 (inverse) decays which are allowed by the thermal masses. 
We consider new production mechanisms which are obtained by including additional 
soft gauge interactions with the plasma. We show that an arbitrary number of such 
interactions gives leading order contributions, and we sum all of them. The rate turns 
out to be smooth in the region where the 1^2 processes are kinematically forbidden. 
At higher temperature it is enhanced by a factor 3 compared to the 1 <H> 2 rate. 
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1 Introduction and motivation 



One of the outstanding problems of standard cosmology is to explain the origin of 
the asymmetry between matter and antimatter. Without such an asymmetry, all the 
structures we observe today would have never formed and mankind would not exist. 
The asymmetry can be expressed as the baryon-to-photon ratio 

— = (6.21 ±0.16) ■ KT 10 (1) 

n 7 

whose numerical value is obtained from a combined analysis of data for large-scale 
structure and the spectrum of the Cosmic Microwave Background [T]. 

In order to obtain a net baryon asymmetry, only three conditions need to be met, 
as outlined by Sakharov in his seminal paper [2]. Yet, providing a model that can 
successfully explain the measured baryon-to-photon ratio remains a challenging task. 
Several different scenarios how to realize the Sakharov conditions have been devised [3] • 
In the last decade, leptogenesis [I] has become very popular. The basic idea of most 
leptogenesis models is to enlarge the particle content of the Standard Model (SM) 
with heavy Majorana neutrinos. In the simplest realization they interact with the SM 
particles via a Yukawa coupling to ordinary, left-handed leptons and the Higgs bosons 
as follows: 

^int = hijNifthj + h.c. , (2) 

where N stands for the Majorana neutrinos, £l and f are the left-handed lepton doublet 
and the Higgs doublet, and cp = Sip* with e a p = —ep ai £12 = 1- Finally, the indices i,j 
label the fermion families, and hij is the Yukawa coupling matrix which need not be 
diagonal. 

The Majorana neutrinos are unstable and decay both into leptons and antileptons, 
N — > £tp, N — > £(p\ The CP symmetry is violated and the corresponding decay 
rates are not equal. Therefore, an excess of antileptons over leptons can be generated. 
The resulting asymmetry is converted into an excess of baryons over antibaryons via 
the sphaleron transitions which conserve B — L but violate B + L [5] . In addition to 
providing a source for the measured baryon asymmetry, this scenario offers a framework 
to explain the smallness of the neutrino masses via the seesaw mechanism [BJ. This 
twofold virtue is what makes the scenario of leptogenesis particularly appealing. 

Despite a substantial amount of work and progress El [91 [10] a complete theory 
of leptogenesis is still lacking [11]. In this paper we study a type of processes which 
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so far has not been considered in this context. We show that they contribute to the 
production of heavy Majorana neutrinos at leading order in the coupling constants, and 
we find that their contribution is numerically large. Our results therefore constitute 
an important step towards a complete treatment of thermal leptogenesis. 

We compute the production rate of a Majorana neutrino in a hot electroweak plasma 
that is fully equilibrated, except for the Majorana neutrinos themselves. The produc- 
tion rate is part of the network of Boltzmann equations that is solved to obtain the 
baryon asymmetry. Since it sets the initial conditions for leptogenesis, it is also of 
practical interest to study the production rate by itself. We assume that the number 
density of Majorana neutrinos is small compared to the equilibrium density, so that 
the inverse processes which reduce their number density can be neglected. We focus on 
the lightest Majorana neutrinos N\ = N, which we assume to be the dominant source 
of lepton asymmetry. 

When the temperature T is sufficiently above the Majorana neutrino mass Mn, a 
peculiar type of production mechanism occurs. It was already considered in Refs. [SUH]- 
Interactions with the hot plasma generate thermal masses, which are much bigger for 
SM particles than for the Majorana neutrinos. Therefore the SM particles can become 
"heavier" than the Majorana neutrino, and the decay of a Higgs boson into Majorana 
neutrino and SM lepton can become possible. Since thermal masses are parametrically 
small compared to the typical particle momentum, all momenta involved in this decay 
are nearly collinear. In this paper we show that there are additional nearly collinear 
processes, involving soft electroweak gauge interactions, which contribute to the leading 
order production rate. 

We focus on the leading order in the SU(2) and U(l) gauge couplings g and g', 
the top quark Yukawa coupling constant h t and the Higgs self-coupling A. We do not 
consider the production via 2 -h- 2 scattering processes Q For the power counting we 
assume that all these couplings are of the same order and collectively refer to them by 
g. All other SM couplings are neglected. We perform the computation in the high- 
temperature regime where <C T. This allows us to formally treat the mass of the 
Majorana neutrino as being soft, ~ gT, and therefore parametrically of the same 
order as the thermal Higgs and lepton masses. We demonstrate that even at leading 
order the production cannot simply be understood in terms of scattering processes 

3 In the literature one can find several calculations of 2 f> 2 scattering rates (see, e.g. [8j [12j [13l 
HUES). However, to the best of our knowledge, there is no calculation which consistently treats all 
leading order thermal effects. 
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involving only a handful of particles. W-strahlung' and inverse decay processes involv- 
ing multiple interactions mediated by soft electroweak gauge bosons are not suppressed 
despite the large number of vertices. The emission occurs almost collinearly, so that 
propagators are nearly on-shell and compensate the suppression. In position space the 
radiated particle and its source overlap over large distances, and the interference of 
different interactions cannot be neglected [16]. This phenomenon has been studied in 
various contexts such as parton energy loss [IT], photon [18], [191 120] and gluon [21] 
production in a quark-gluon plasma (for a general discussion see [22]). Recently [2U] . 
we presented a new approach how to consistently include soft gauge interactions in the 
computation of a thermal particle production rate. It is formulated in a way that is 
largely independent of the type of particles whose production we want to study, and 
can therefore easily be adapted to the case at hand. 

The paper is organized as follows. In Sec. [2] we relate the production rate, which 
describes out-of-equilibrium physics, to a real-time correlation function in equilibrium. 
The latter can be calculated in thermal field theory, which is subject to the following 
sections. We describe the physics of collinear emission and outline the relevant momen- 
tum scales in Sects. 13.11 and 13.21 In Sec. 13.31 we give a short and qualitative summary 
how we proceed to obtain the leading order production rate due to collinear emission 
processes. The rest of Sec. [3] provides all the technical details that are needed to arrive 
at the final results, and the reader who is not interested in the details of their derivation 
may skip directly to Eq. (139]) . In Sec. H]we present numerical results and we conclude 
in Sec. [5j The appendix finally explains how to obtain the numerical solutions. 

2 Production rate and thermal field theory 

At lowest order in their Yukawa couplings the production rate of Majorana neutrinos 
with 4-momentum k can be written as n 




Here S ret (A;) = S(/c° + i0 + , k) is the retarded self-energy for the Majorana neutrinos, 
and /p is the Fermi-Dirac distribution. 

4 Our formula for the production rate of a spin 1/2-fermion is consistent with those shown in the 
literature (e.g. p!3,) although the overall sign appears to be different. However, with our conventions 
the self-energy corresponds to (—1) times the Feynman diagrams. 
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Figure 1: The imaginary part of this diagram gives the production rate due to decay 
and inverse decay processes. This diagram also determines the inhomogeneous term in 
Eq. (T23D. 

Due to the Majorana-nature of N there are two types of diagrams contributing 
to the self-energy which differ by the orientation of the internal fermion lines. We 
neglect Standard Model CP-violation. Then both types of diagrams give the same 
contribution. Therefore we consider only one of them, and multiply by 2 to obtain the 
correct rate. We do this by considering only the self-energy for right-handed TV's, so 
that the non- vanishing components of the self-energy fit into a 2x2 matrix which we 
denote by £r. Then the rate can be written as 

dr 2 

¥k = j2^ Mk ° )Tl P- kl ™^Mk)] , (4) 
where a = 1, <x = —a. 

3 Production by collinear emission 
3.1 Production processes 

Without any SM interactions the imaginary part of the self-energy is obtained by 
cutting the diagram in Fig. [TJ with the lines in the loop being tree-level propagators. 
The corresponding production mechanism is the inverse decay of the massive Majorana 
neutrino, tp£ — > N. 

The typical momentum of a particle in the hot plasma is 0(T). We refer to these 
momenta as 'hard'. In the high-temperature regime, when <gT, one cannot ne- 
glect the modification of the dispersion relation of SM particles, which is caused by the 
interactions with other particles in the hot plasma. The dispersion relation for hard 
particles can be written as p 2 = m 2 , where m is the so-called asymptotic mass [23]. % 

It is oftentimes referred to clS 777. oo ■ 
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Figure 2: Decay and inverse decay processes that contribute to the production of 
Majorana neutrinos. 

It is given by the self-energy of a hard particle with (nearly) light-like momentum, and 
the dominant contribution is due to hard loop momenta. Its values for the Higgs and 
the lepton doublets are given by§ 



Note that the gauge field contributions to the asymptotic masses for Higgs and leptons 
are equal (cf. Refs. [22]). However, the Higgs also receives important contributions 
from the Yukawa interaction with the top quark and from the Higgs self-interaction, 
so that m v > me. All other contributions can be neglected due to the smallness of 
the corresponding coupling constants. Furthermore, the thermal mass of the Majorana 
neutrinos can be neglected, so that M N has a temperature-independent value. There- 
fore, the Higgs and the leptons become "heavier" than the Majorana neutrinos at high 
T. 

When m y + me > the inverse decay is kinematically forbidden. At even higher 
temperature, when m^ — me > Mn, the phase space for the decay <p — > £N of the Higgs 
boson opens up (see Fig. [2]). The decay due to thermal masses has been considered as 
a mechanism for Majorana neutrino production previously [HI E] - 

When M^r < gT, the characteristic feature of both processes described above is that 
the momenta of the particles are nearly collinear because the relevant masses are small. 

6 For scalars the asymptotic mass is the same as the thermal self-energy computed for vanishing 
momentum, which enters the finite temperature effective potential for the scalar field. It also equals 
the frequency of scalar field oscillations with zero p. For fermions, however, the asymptotic mass is 
larger than the oscillation frequency for vanishing p by a factor \[2 [53] . 

7 In Refs. [H [9] the oscillation frequency for vanishing p was used instead of the asymptotic mass 
for fermions. 





^(3g 2 + g' 2 + 4h 2 + S\) T 2 , 

rA^ + 9 ' 2 )T 2 . 
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Figure 3: Example for two processes whose interference needs to be taken into account 
in a consistent leading order treatment. 

The angles between the particle momenta are of order g. We will show that multiple 
scattering mediated by soft gauge bosons contributes already at leading order, similar 
to the thermal production of on-shell photons [HI [19] . The additional couplings are 
compensated by nearly on-shell propagators. The lepton and the Higgs boson undergo 
an arbitrary number of scatterings off soft gauge bosons during the 'emission' of the 
Majorana neutrino which still feels the presence of its 'source'. In order to take this 
phenomenon into account, one has to resum an infinite set of diagrams, like in [T9| [20] 
for the production of photons from a quark-gluon plasma. This will be dealt with in 
the rest of Sec. [3j 

3.2 Momentum scales 

For the processes that we consider the 4-momenta of the emitting and of the emitted 
particles are almost collinear and close to the light-cone. The 3-momenta point in 
approximately the same direction, which we represent by the 3- vector v with v 2 = 1. 
The components parallel to v are denoted by 



We further define the light-like vector v = (l,v). 

One has to account for three distinct momentum scales: 

1. The emitting particles and the emitted particle both have p\\ ~ T, which is our 
hard scale. 

2. All 3-momenta perpendicular to v are soft, p± ~ gT. Furthermore, all momen- 
tum components of the exchanged gauge bosons are soft, ~ gT. 

3. Finally, all 4-momenta k satisfy v ■ k = (ko — k») ~ g 2 T. 

The hard loop momenta have k 2 ~ g 2 T 2 . Therefore the propagators are sensitive 
to the asymptotic mass m ~ gT. 



P\\ =pv. 



(6) 
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Figure 4: Example for a self-energy diagram whose imaginary part contributes to the 
leading order production rate. 

3.3 Strategy of the calculation 

We now turn to computing the production rate due the processes discussed in Sec. 13.11 
They involve interference terms with an arbitrary number of soft SU(2) and U(l) gauge 
bosons and the corresponding self-energy needed for ([3]) is given by the sum of all ladder 
diagrams of the type shown in Fig. HI We use the approach of Ref. [20] to obtain an 
integral equation for the self-energy that we can subsequently solve numerically and 
get the contribution to the production rate. Before going into the details, let us recall 
the strategy of the calculation [20] : 

1. Integrate out the hard field modes. This generates the asymptotic masses for 
hard particles near the light-cone, i.e. with p 2 ~ g 2 T 2 . Furthermore, one obtains 
the Hard Thermal Loop (HTL) effective action for soft fields, in particular the 
gauge fields. At leading order no thermal width for the hard particles is generated 
in this step. The width due to hard-hard interaction is of order g A T and can be 
neglected. The leading order width is due to soft gauge interactions and enters 
only at a later stage. 

2. Consider 1-loop diagrams with 2 external Majorana neutrinos and with an arbi- 
trary number of external soft gauge bosons, in the limit that the loop momentum 
and the hard external momenta are nearly collinear. Expand the propagators and 
vertices in the coupling and keep only the leading order terms. First compute 
the 2-point function without any gauge boson explicitly. Then set up a recursion 
relation that relates a given n-point function to (n — l)-point functions where 
one of the gauge bosons has been removed. This recursion relation is most easi- 
ly formulated in terms of a 'current' which is induced by background Majorana 
neutrino and gauge fields. 

3. Finally, integrate out the soft gauge boson background. The gauge bosons then 
only appear in self-energy insertions (thereby generating the thermal width for 
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the lepton and the Higgs boson) and as rungs in the ladder diagrams (see Fig.HJ.Lj 
The new current satisfies an integral equation that is straightforwardly obtained 
from the one found in step 2. By stripping off the external field we obtain an 
integral equation for the desired self-energy 

In the remainder of Sec. [31 we will provide all the relevant details to arrive at the 
final result fl39l) for the production rate of Majorana neutrinos, following the strategy 
outlined above. 



3.4 2-point function 

In this section we explicitly compute the 2-point function of the Majorana neutrinos 
shown in Fig. [TJ ignoring any interactions, except that we use the asymptotic masses. 
It turns out that all n-point functions with additional {n — 2) soft gauge bosons can 
be obtained from (n — l)-point functions through a simple recursion relation. All 
simplifications of propagators and vertices which are needed for the general case already 
appear in the calculation of the 2-point function. 

Let us first have a look at the propagators. Keeping only the leading order g 2 T 2 
terms in the denominator, the scalar propagator with mass m v can be written as [2U] 

A (P) = ~2 2 - O 7 

p 2 - ml 2p|| 

with 

DaijP) = , 2 ~]_ 2 , /(0 y (8) 

Now consider a spin 1/2 propagator. We consider only left-handed fermions prop- 
agating in the loop. Thus we can work with 2-component Weyl spinors. One can 
approximate 

S(p)^^a-p (9) 
2 P\\ 



with 



p = p--^u, (10) 
2p|| 



8 At leading order, only diagrams of the form shown in Fig. 2] are generated. No diagrams with 
crossed ladder rungs or nested loops occur |19l 120] . 
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where u = (1,0) is the 4- velocity of the plasma. The matrix a ■ p contains terms of 
order T, gT and g 2 T\ higher order terms have been neglected. It is necessary to keep 
the g 2 T terms because they give a leading order contribution to the rate in Eq. ([3]). 

If the fermion momentum is on-shell, p 2 = mj, then p is light-like up to higher 
order terms which we have neglected. Therefore, if we evaluate the loop integral in 
the imaginary-time formalism and take the imaginary part, the fermion propagator is 
on-shell, and we can treat p as a light-like 4-vector. This allows us to write 

a-p = 2p\\r){p)ri\p), (11) 

with a 2-component Weyl-spinor r)(p). Then the fermion propagator becomes 

S(p)=r ] (p)rf(p)D e (p). (12) 

We choose the 3-axis in the direction of v and expand r](p) in powers of g, 

V = Vo + Vi + 0(g 2 ) (13) 

with 

*>=(?). m = -^po. (14) 

The lower component is 0(1) and the upper component is 0(g). Since we only need 
the leading order for each component of r), which is at most of order g, we can write 
rj (p) instead of r\ (p) because the difference of p and p is of order g 2 T. 
We use the partial fractioning 

D e (p)D^(p - k) = — — [D e (p) - D^p - k)] , (15) 
e(k,p) 

where 



(P± ~ k ± ) 2 + ml pi 



me 



e(t,p) 5 ».t + - (p - (16) 

is the difference of the energy poles of the two propagators. Then it is straightforward 
to compute the self-energy. It will turn out to be convenient to write it in the form 

Z R (k) = \h\ 2 J ^r,(p)S(k,p). (17) 
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Figure 5: 1-loop diagrams with soft external gauge field lines. Only the 2-point function 
Fig. [1] needs to be calculated explicitly. The n-point functions with n > 2 are related 
to the (n — l)-point functions by a recursion relation. 



It is proportional to \h\ 2 = Ylj \h\j\ 2 where the sum goes over the families of leptons. 
The gauge interactions which will be included are the same for each lepton family, 
so that this factor can always be extracted. Our result for the 'reduced self-energy' 
T,(k,p) is 

2e(k,p) P\\-k\\ 

where d(r) = 2 denotes the dimension of the gauge group representation for lepton and 
Higgs, and 

&(k h p\\) = fj,(jp\{) + f B (p\\-k\\). (19) 

From Eq. (Tl8~I) one then obtains the rate due to the decay of the Higgs or from the 
inverse decay of the Majorana neutrino, with the diagrams shown in Fig. |2j Explicit 
results can be found in Sec. 13.81 



3.5 Self-energy in soft gauge field background 

Now we consider diagrams with 2 external iV's and, in addition, (n — 2) external gauge 
field lines. These diagrams can be recursively related to diagrams with {n — 3) gauge 
fields. To obtain this relation we pick out one vertex, the leftmost one in Fig. [51 with 
incoming momentum k. As for the 2-point function we can use the approximations 
(jSJ) and f ll2p for the propagators. After the partial fractioning ( |15p this diagram is 
proportional to the difference of the diagrams in which either the external line with 
momentum q or the one with momentum q' is omitted. 

It is convenient to work with the generating functional of the diagrams rather than 
with the diagrams themselves. We attach background fields to the external lines, 
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except to the one on the left. That means that we consider the first derivative of the 
generating functional. It can be interpreted as a 'current' which is generated by the 
presence of the background fields. 

Here we simplify the treatment of Ref. [20], where first diagrams were considered 
and then the generating functional was built from them. Instead we start directly 
with the current. We define J(k,p) as the sum of all diagrams in the background of 
external Majorana neutrino N and of gauge fields A£, but without the integral over 
the 3-momentum p, without the trace over the SU(2) indices, and without the vertex 
factor rj(p) (cf. Eq. (j!7j) ). Each lepton family gives the same contribution. Therefore 
we can also pull out the factor \h\ 2 as in Eq. (jl7|) . The current without gauge fields, 
which corresponds to our result from Sec. 13.41 is 

tTJ(k,p) = E{k,p)N(k). (20) 

Here 'tr' refers to the trace over SU(2) indices. 

For computing higher n-point functions it is furthermore convenient to associate 
the factor (2p||) _1 in the scalar propagator A(p) (see Eq. (J7j)) with the vertex in the 
direction of p, rather than with the propagator itself. Then the vertex of two Higgs 
lines with one gauge field with soft momentum q is 

-L(2p-qr = v» + 0(g). (21) 
2 P\\ 

On the fermion lines we associate the spinors r)(p) and rf(jp) not with the propagator, 
but with the vertices to which the propagator is attached. For a fermion-gauge field 
vertex one therefore obtains a factor 

v\p-q)^v(p) = v" + 0(g). (22) 

From Eqs. ()2ip and (T22]) we only need the leading order contributions, which are equal. 

For the two propagators which are connected to the left vertex in Fig. [5] we use the 
same approximations (jHJ) and ( 112]) as for the computation of the 2-point function, and 
we use the partial fractioning ([15]) . The resulting two terms are then proportional to 
diagrams in which either the propagator with momentum p or the one with momentum 
p — k has been omitted from Fig. Thus, if one also leaves out the vertex factors 
and the background gauge fields connected to these propagators, each of the two terms 
gives a (n— l)-point function. For the second term in Eq. ffl5]) we obtain a contribution 
in which the propagator with momentum p has been omitted. We perform a shift in 
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the summation variable, p° — > p° + q°- Then the remaining propagators are the same 
which appear in the (n — l)-point function, but with the loop 3-momentum p replaced 
hy p — q. Therefore we can write J(/c, p) in terms of the difference of J times a gauge 
field and of a gauge field times J, 

e(k,p)J(k,p) = - \^{k h p^^-N{k) 

Z T)\\ K\\ 



+ 



J(k-q,p)vA(q)-vA(q)J(k-q,p-q) . (23) 



In the imaginary-time formalism 

g u =tui 

where a; is a bosonic Matsubara frequency. Furthermore, A^ = A^T A contains both 
SU(2) and U(l) gauge fields. Note that we have included the gauge couplings in A A . 
They will appear in the gauge field propagators. Eq. (|2"3"|) is of the same form as Eq. (21) 
in [20], only the inhomogeneous term on the RHS is different. 

3.6 Integrating out the soft gauge bosons 



Now we use Eq. (123]) to integrate out the soft gauge fields. This will give us the N 
self-energy including soft gauge interactions. To see how it works, we write Eq. (|23|) 
schematically as J = N + AJ, leaving out all terms and all factors which are not 
essential for this purpose. We iterate this equation once, which gives J = N + A(N + 
AJ). Now we integrate this expression over the gauge fields. The term linear in A drops 
out and we obtain / = N + (^AAJ^, where (■ • • ) denotes the path integral over the 

soft A-field. In [20] we have shown that at leading order (^AAJ^ = i^^j m ^ s 



way we obtain a closed equation for ^ J ^ 

Now we can become more explicit. We iterate Eq. (|23|) once and we take the trace 
over the SU(2) indices. This allows us to bring all terms in the same order, 



e(k,p) tr (J(k,p)) = -^&(k hPl {)^^N(k) 



1 



v ■ (k — q) 



tr 



vA{q)vA{q'))(j{k,p)- J(k,p hP± -q x )) . (25) 
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Inside the integral we have approximated e(k — q,p) ~ e(k — q,p — q) ~ v ■ (k — q), 
i.e. we have neglected the terms containing transverse momenta and thermal masses, 
even though they are of the same order as the term we kept. This is possible because 
the gauge field propagator separately depends on q ~ gT and q\\ ~ gT, and not on 
their difference v ■ q ~ g 2 T. Therefore the terms we have omitted only contribute to a 
higher order shift of the integration variable q» . 

For the SU(2)xU(l) gauge fields one has at leading order 

(MdMfl) = ha + q') [CiWA^q) + y y 2 ^M] (26) 

where in the imaginary-time formalism 

S( q + q>) = T-% 0+q , j0 (27r) 3 S(q + q'). (27) 

Furthermore, C^r) is the quadratic Casimir operator which for the fundamental re- 
presentation of SU(2) equals 3/4, and y% = —1/2 is the lepton hypercharge. Finally, 
A^ u and A'^ are the HTL resummed propagators of the SU(2) and U(l) gauge fields. 

On the RHS of Eq. (|25|) the unknown function J does not depend on q or q\\. 
Therefore these integrals can be performed. One encounters the expression 



and similarly with A' . Here k is purely imaginary, k = iu, where Co is a fermionic 
Matsubara frequency. After performing the sum over the bosonic Matsubara frequen- 
cies u, one has to analytically continue to k° + i0 + , where k° is now real. Using 
standard results for the HTL resummed propagators [26] and the sum rule of Ref. [27] 
one obtains 

I(k° ± z0 + , k, qx) * T-T (\ - 2 ] 2 ) (29) 

2 W± Ql + m bJ 

where ttid is a Debye mass. In the SM the Debye masses of the SU(2) and U(l) gauge 
bosons m D and m' D are given by [2S] 

2 11 2rr2 l 2 11 /2rr2 / on \ 

m D = Y 9 T ' m D = Y 9 T ■ ( 30 ) 

The RHS of Eq. ( 129]) no longer depends on k. The only /c-dependence of I is through 
the sign of the imaginary part of the frequency. 
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Now we strip off the background iV- field from J which gives £ (see Eq. fl2"Uj) ). The 
equation for T,(k,p) then reads 



ie(k,p)Z(k,p) 



2 p\i — K\\ 



+ 



dV 



tf(q ± ) E(k,p)-?:(k,p h p ± -q ± ) 



with the kernel 



C 2 (r)g 2 ( ~^ 



vh' 2 \ -2 



1 



i 



/ 2 



Q± Q± + m D 



(31) 



(32) 



The integral equation ( l3Tj) is of the same form as the one for the production of 
photons (20], only the inhomogeneous term is different. Both e(k,p) ~ g 2 T and 
d 2 qj_ ( if(q±) ~ g 2 T. Therefore the integral term which resums the multiple interac- 
tions with an arbitrary number of soft gauge bosons is of the same order as the term 
which corresponds to the tree-level (inverse) decay. 

3.7 Production rate of Majorana neutrinos 

To compute the trace over spinor indices in Eq. @ we expand 



a-k = 2/qiXoxJ) 
where t]q is given by Eq. (fi4|) . and 



M 2 



N Vovl + 0(g i T) 



2h 



(33) 



Xo 



(34) 



Since the first term in Eq. (|33|) is 0(T) and the second is 0(g 2 T), both terms in 
Eq. (Tl3|) . and consequently both components of £(/c,p), contribute to the leading 
order production rate. 

Now, if one has the solutions to the equations 



H k ,P)f(P±) 
ie(h,p)if>(p x ) 



dV 

(27T) 2 

dV 
(2vr) 2 



V(q±) if(p±) ~ f(P± ~ q±)] = 2pj 



%i)Wpi)-^(Pi-9i)] 



(35) 
(36) 
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which were obtained in [T^l ED] for transverse and (virtual) longitudinal photon pro- 
duction, then the solution to Eq. (13T1) is given by 



E(fc,p) 



2 P\\ - 



>P\\ 



(37) 



Unless one ignores the interactions with soft gauge fields by putting ^ equal to zero, 
the imaginary part of k° is no longer relevant in Eqs. ( 1315]) and (1361) . Only its sign has 
entered by determining the sign of c & through Eq. (1291) . 

We may now choose v in the direction of k, so that = 0, and kn = |fe| = k. 
Combining Eqs. Q, flUD, (ffB> O, (E7J) and using 



/f(*)[/f(p) + HP ~ k)] = -/ F (p)/ B (fc - p) 



(3? 



we obtain for the production rate 

dr c/(r)|/i| 2 /• d 3 p 1 



d 3 fc 



(2tt) 3 2A; J (2tt) 3 A; — p. 



/f(P||)/b(^ — P||) Re 



k n M 2 N 

7T^P± • / + -r-4> 



2p 



(39) 



The integral equations ( )35|) and (136|) for / and ^ can only be solved numerically 
which is greatly simplified by transforming them via a Fourier transformation into a 
boundary- value problem for a differential equation as shown in appendix [A] 

The Yukawa interaction couples Majorana neutrinos to SM leptons with opposite 
chirality. In the massless limit where helicity equals chirality the collinear emission 
is forbidden by angular momentum conservation. This argument is not affected by 
the thermal masses for fermions because they do not violate chiral symmetry. The 
contribution containing / corresponds to helicity changing processes, and therefore it 
vanishes when p± is zero, x/i does not vanish in the collinear limit and corresponds to 
helicity conserving processes. Therefore ip vanishes in the limit — > 0. 



3.8 Tree-level (inverse) decay 

In addition to solving the full integral equations fl35|) and fl36|) . we want to explicitly 
calculate the contribution from the tree-level processes shown in Fig. [2]in order to study 
how strongly the multiple soft scattering affects the production rate. This corresponds 
to neglecting the integral term in (I3T1) completely, and the equation becomes a purely 
algebraic one. In this case one has to keep a small imaginary part of k° in e(k,p). 
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Taking the real part of the trace in (j2HD and performing the integration over p± yields 



p+ 

dr tree ^ d(r)\h\ 2 f dp\\ 



■/p(P||)/b(* -P||) [p\\{M%-ml) - (k-p {l )ra 2 ] . (40) 



d 3 k 4(2tt) 3 A; 2 J 2n \p 
p- 

The integration limits are determined by the conditions that the energy of the 
lepton be positive and that the particles be on-shell, which means that e(k,p) = has 
a real solution for \p± | . For the two processes shown in Fig. [2J this gives the integration 
limits 

x± Vy 

where 



P ± = ^2MT (41) 



X = M 2 N + ra 2 - m%, 

Y ee (rritp + m l + M^lm^ -m £ + Mjv)(m v + m e - M N )(m v - m e - M N ). (42) 

It is easy to see that Y > only if either m v > M^r + (which corresponds to 

the decay of the Higgs boson) or Mn > rn^ + rri£ (which corresponds to the inverse 

decay of the Majorana neutrino). For — me < Mn < rn^ + ran we therefore obtain 
dr trcc / / d 3 A , = Q 



4 Numerical results 

To obtain numerical results we have to specify several parameters. The mass Mjq 
of the Majorana neutrino and its Yukawa couplings are unconstrained by low-energy 
neutrino physics. We have chosen the exemplary value Mm = 10 7 GeV and we always 
plot our rates divided by \h\ 2 = J2j \^ij\ 2 - The SM couplings are evaluated at the 
scale /i = 2ttT using the 1-loop renormalization group equations [29]. To determine 
the Higgs self-coupling A we have assumed ran — 150 GeV for the zero temperature 
Higgs mass. 

First consider the integrated production rate T = J d 3 k(dT /d 3 k) as a function of 
z = Mn/T. Our results are shown in Fig. [6j At small z the 'tree-level rate' (I40p is 
due to the decay of Higgs bosons, while at large z it results from the inverse decay of 
heavy Majorana neutrinos. At intermediate z none of these processes is kinematically 
allowed, and the rate vanishes. It is remarkable that the full line which includes soft 
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Figure 6: Number of produced Majorana neutrinos per unit time and unit volume 
as a function of z = Mjq/T. The dotted curve is the result without any soft gauge 
interactions. The full line includes an arbitrary number of soft gauge interactions. 

gauge interactions is very smooth in the regions where the decay of the Higgs boson 
becomes kinematically forbidden. 

One can see that the full rate is larger than the tree-level rate by about a factor 3 at 
small z. It decreases only mildly when the tree-level processes are forbidden. When the 
inverse decay process sets in, the difference between tree-level rate and the complete 
rate goes to zero. This is expected, since the collinear enhancement is a relativistic 
effect and it disappears when the Majorana neutrinos become non-relativistic. One 
should emphasize that the strong enhancement caused by the soft gauge interactions 
does not signal a breakdown of perturbation theory, because all contributions discussed 
above are leading order. 

It is also interesting to consider the contribution due to helicity changing and helicity 
conserving processes separately (cf. the discussion at the end of Sec. I3.7|) . The rate of 
helicity changing processes does not vanish in the limit Mjv — > 0, and should therefore 
be dominant at small z. The results are shown in Fig. [7J We clearly see that the 
helicity changing process dominates at high temperatures and the helicity conserving 
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0,01 



helicity-flip, tree-level 

— helicity-flip, full 

— non-flip, tree-level 
non-flip, full 




z = M N /T 

Figure 7: Contributions to the production rate of Majorana neutrinos due to helicity 
changing and helicity conserving processes as a function of z = Mn/T. 

process is relevant only at low temperatures, T ~ Mat, while it becomes negligible at 
T Mat. This is true both for the tree-level processes and the processes which include 
multiple soft scattering. 

Finally we study the momentum spectrum of the produced heavy neutrinos. If it 
was thermal, the rate would be proportional to the Fermi-Dirac distribution /f(&). In 
Fig. [S]we show the ratio of the differential production rate and /f(&) for two different 
temperatures. For z = 0.1 the production occurs both via a tree-level process (in 
this case, the decay of the Higgs boson) and via processes involving multiple soft 
scattering whereas for z = 0.5 the tree-level processes are kinematically forbidden. We 
see that the spectrum is peaked in the infrared which means that the typical momenta 
of the produced neutrinos are smaller than in equilibrium. This holds both for helicity 
changing and for helicity conserving processes. 

5 Summary and Conclusions 

In this paper we have studied the production of relativistic Majorana neutrinos, rel- 
evant for models of thermal leptogenesis, in a hot electroweak plasma. Based on our 
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Figure 8: Differential production rate divided by Fermi-Dirac distribution function for 
two different temperatures. 

previous work [20], we have obtained an equation which sums all leading order collinear 
production processes. It includes the tree-level processes-the decay of the Higgs boson 
as well as the inverse decay of the Majorana neutrino-and in addition processes involv- 
ing multiple scattering mediated by soft gauge boson exchange. All these turn out to 
contribute at leading order in the coupling constants and have never been included in 
previous treatments of thermal leptogenesis. 

Numerically, we find a very pronounced increase in the thermal production rate 
when soft gauge interactions are included. At high temperatures, when the production 
via tree-level Higgs decay is allowed, the rate increases by about a factor 3. When the 
tree-level processes are forbidden, the rate in units of temperature drops only slightly, 
and the production of Majorana neutrinos remains effective. When the temperature is 
close to Mn the rate is not significantly affected by soft gauge interactions. We further 
showed that the production rate is dominated by helicity changing processes and that 
helicity conserving (inverse) decays only play a significant role at low temperatures, 
T ~ Mjv- In addition to the production rate, we have studied the momentum spectrum 
of the produced Majorana neutrinos and found that it is strongly peaked in the infrared. 
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The spectrum is thus not even approximately thermal. 

The leading order production rate of heavy Majorana neutrinos also receives con- 
tributions from 2 <H- 2 scattering processes. They have been considered previously in 
the context of thermal leptogenesis (see e.g. [8j [121 H3J HH US]), but a complete and 
consistent leading order calculation has not been done so far. In [8] the Higgs decay 
was found to dominate over the 2 H 2 processes at high temperature. In this paper 
we have shown that adding soft gauge interactions to the Higgs decay process leads to 
a strong enhancement of the production rate. Therefore it would be very interesting 
to see how the processes discussed here, together with the 2 <H- 2 scattering processes 
affect various scenarios of leptogenesis. 

Acknowledgments This work was supported in part through the DFG funded 
Graduate School GRK 881. 

A Solving the integral equations 

The p_i_-integrals of the functions / and ip in Eq. fl39|) can be obtained without solving 
the integral equations (T53]) and This is achieved by Fourier transformation which 
turns Eq. (|5o^) and (1361) into boundary value problems. Our approach which is described 
in this Appendix closely follows Refs. [27J ED]- 

A.l From integral equation to differential equation 

Eqs. fl35|) and fl36l) are of the same form as the integral equation for the production of 
transverse and of longitudinal photons. Following [27J [3U] we Fourier transform § 



Then the integral over the perpendicular components in (I39I) can be replaced by the 
limit of their Fourier transforms for b — > 0: 




(A.l) 




(A.2) 



(A.3) 



9 We use the same symbol for the functions /, %j) and their Fourier transforms. 
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The integral equations (|35p . f[3T)j) now turn into differential equations for f(b) and ip(b). 
With k± = and k 2 = M 2 N we can write 



e(k,p) = p(p 2 ± + M 2 s ) 



(A.4) 



where 



A' 1 1 



_ M 2 _ PjKgii - k \\) M N - k \\(n - k \\) m e + hPfK 1 ^ A5 ^ 

2P||(P|| - fc||)' eS kfi 



Then we obtain differential equations which closely resemble those in [30 



where b = \b\, and 



with [23 ED] 



(A - M e 2 ff ) /(ft) - JT(&)/(6) = - i2V5(b) 
-H3 (A - M e 2 ff ) ^(6) - JT(6)^(6) = 5(b) 



X(b) = T C 2 (r)g 2 D(m D b) + y 2 g' 2 D (m' B b) 



D(x) 



1 

2^ 



x 



7e + In - + K (x) 



(A.6) 
(A.7) 



(A.8) 



(A.9) 



Here, 7e = 0.5772... denotes the Euler-Mascheroni constant and Kq(x) is a modified 
Bessel function. 

The convergence of the Fourier integrals enforces the boundary conditions 



lim /(b) = 0, lim ip(b) = 0. 



(A.10) 



In the limit b — > all terms which contain no derivatives can be neglected in (1A.6[) . 
(1A.7D . and we obtain the limiting behaviour 



f(b)=c f - + 0(b), 



(A.11) 

4j{b) = c^ In b + 0(6°). (A.12) 

The constants Cf, are fixed by the 5-functions in Eqs. ( IA.6j) . ( 1A.7|) which gives 

1 



c 4' 



71(3 ' 

i 

2t7p 



(A.13) 
(A.14) 
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Due to rotational invariance we must have f(b) = h(b)b and ip(b) = ip(b). In terms 
of the new function h(b) the RHS of ( 1 A. 21) becomes 2 lim^o im h(b). We numerically 
solve the ODEs for h(b) and ip(b) 

-i/3 (dl + ^d b - h(b) - X{b)h{b) = (A.15) 

-i/3 (d 2 b + ^d b - M 2 c ^j fab) - X{b)fab) = (A.16) 

for b > 0. The conditions (1A.10p . ( 1A.13|) and (1A.14|) are then sufficient to determine 
the solutions unambiguously. 

A. 2 Numerical procedure 

Finally we want to describe our algorithm that we used to obtain the solutions shown 
in Sec. HI Both equations (IA.16I) and (IA.15[) can be solved by the same method and 
we only describe the procedure in terms of flA. 16h explicitly. 

For the numerical solution, it proves convenient to split the function into a tree-level 
part and another part coming from multiple soft scattering, 

fab) = fa(b) + fa(b). (A.17) 

The function ipo solves flA.7j) with D(b) = and can be written in terms of Bessel 



functions. For b — > 0, the function D(b) behaves like D(b) ~ b 2 lnb. This implies that 
the limiting behaviour of i/jq for & — ?- is also given by Eq. (1A.14D . Therefore ipi(b) 
must be regular at b = 0. The general solution for fa can be written as 

faib) = Cl ^\b) + c 2 ^\b) + ^(b) (A.18) 

where ip^^ is a particular solution of the inhomogeneous equation while the linearly 
independent solutions ip^ solve the corresponding homogeneous equation. One can 



choose ip^ such that it has the limiting behavior of Eq. (1A. 12j) . and a ip^ which is 



regular for b — > 0. Thus we must have c\ = 0. The other constant is then found by 
imposing flA.101) and we obtain 

c 2 = - lim C^). (A.19) 

Numerically, one has to choose finite values for '6 — > 0' and '6 — > oo\ We have found 
bo = 10~ 5 T~ 1 and b^ = 30T -1 to be a reasonable choice, giving stable numerical 
solutions. 

The algorithm providing us with the desired solution is then the following: 
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1. Solve the homogeneous ODE for tpiip) with ipi(b ) = l,ip[(b ) = to obtain 
^ 2) (6), 

2. Solve the inhomogeneous ODE for ipi(b) with V'i(^o) = ^V'i(^o) = to obtain 
^i part) (&), 



3. Compute C2 via (1A.19j) . Then, due to the initial conditions chosen in the previous 



steps, the desired result is given by Re^i(6 — > 0) = Rec2. 

In order to find Im hi(b —> 0) we can use the same algorithm, only in the second step it 
is more convenient to choose /ii(feo) — 1> h[(bo) = because then Imhi(b — > 0) = Imc2. 
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